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Summary. — The paper discusses dynamics of quantum measurements in meso- 
scopic solid-state systems. The aim is to show how the general ideas of the quantum 
measurement theory play out in the realistic models of actual mesoscopic detectors. 
The two general models of ballistic and tunneling detectors are described and studied 
quantitatively. Simple transformation cycle demonstrating wavefunction reduction 
in a mesosocpic qubit is suggested. 



1. Introduction 



Despite the long history of quantum mechanics, the quantum measurement problem 
continues to attract interest motivated mostly by the counter-intuitive features of the 
"wavefunction reduction" . Although dynamics of any measurement set-up is governed 
by the Schrodinger equation with an appropriate Hamiltonian, full description of the 
measurement process can not be obtained without account of the changes in the wave- 
function of the measured system caused by the random process of selection of one specific 
outcome of measurement out of the range of possible outcomes. This selection process is 
trivial in the case of classical dynamics, when all possible outcomes of measurement are 
"orthogonal" and the observation of the measured system in one particular state does 
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not imply any changes in the system beyond the statement that it occupies this and 
not any other state. For a quantum system, however, existence of the non-commuting 
observables implies that selection of one particular outcome of measurement can change 
the wavefunction of the system in a highly non-trivial way. Such a reduction of the 
wavefunction appears as an evolution principle additional to the Schrodinger equation. 
Moreover, the changes in the wavefunction of the measured system induced by it can 
violate the basic features of dynamics which follows from the Schrodinger equation, de- 
spite the fact that the measurement process as a whole is governed by this equation. The 
best known example of this situation is the case of EPR correlations 1 between the two 
spatially separated spins, which violate the no-action-at-a-distance principle as quanti- 
fied by the Bell's inequalities 2. From the perspective of the wavefunction reduction, 
the EPR correlations appear as a result of selection of one random specific outcome of 
the local spin measurement. On average, there is no action-at-a-distance in a sense that 
the correlations by themselves can not lead to information transfer between the points 
where the spins are located. 

Current interest to the solid-state quantum information processing (see, e.g., the re- 
views EI El) j motivates development of mesoscopic solid-state structures that can serve 
both as simple quantum systems, e.g., qubits or harmonic oscillators, and the detectors. 
Although in experiments, the mesoscopic detectors did not reach the stage yet where they 
can be used to look into the basic questions of the quantum measurement theory (which 
requires the quantum-limited detection) one can expect this to happen quite soon. A new 
element introduced by the mesoscopic structures in the discussion of quantum measure- 
ments is the fact that the wavefunction reduction is not necessarily caused by interaction 
of a "microscopic" measured system with the "macroscopic" detector. In mesoscopic 
structures, the measured systems and the detectors are similar in many respects (includ- 
ing dimensions and typical dynamics) and are frequently interchangeable: a measured 
system in one context can act as the detector in the other, and vice versa. This shows 
that the boundary at which the quantum coherent dynamics should be complemented 
with the wavefunction reduction is not universal. 

The aim of this work is to provide a quantitative discussion of models and measure- 
ment dynamics of the mesoscopic detectors. The discussion emphasizes the interplay 
between the dynamic and information sides of the measurement process and can serve 
as an introduction to the problem of wavefunction reduction in mesoscopic structures. 

2. Measurements dynamics of ballistic mesoscopic detectors 

Majority of the mesoscopic detectors use as their operating principle ability of a mea- 
sured system to control the transport of some particles between the two reservoirs. The 
information about the state of the system is contained then in the magnitude of the 
particle current between the reservoirs which serves as the detector output. In the most 
direct form, this principle is implemented in the quantum point contact (QPC) detec- 
tor [SI El) which presently is the main detector used for measurements of the quantum 
dot qubits El HU H3 • In the QPC detector (Fig. the propagating particles are 
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Fig. 1. - Schematic of the QPC detector measuring charge qubit. The two qubit states \j), 
j = 1,2 are localized on the opposite sides of a tunnel barrier and are coherently coupled by 
tunneling across this barrier with coupling strength A. Transfer of the qubit charge between 
the states \j) changes electrostatic potential in the scattering region of the QPC affecting the 
current I through it that is driven by the applied voltage V. 



electrons which move ballistically through a short one-dimensional constriction formed 
between the two electrodes of the QPC. The electrodes can be viewed as reservoirs 
of independent and effectively non-interacting electrons. The measured system creates 
electrostatic potential that makes the scattering potential Uj (x) for electrons in the con- 
striction dependent on the state \j) of the system, and in this way controls transmission 
probability of the QPC. The output of the QPC detector is the electric current I driven 
by the voltage difference V between the electrodes. The current depends on the electron 
transmission probability, and as a result, contains information about the state \ j). Since 
interaction between the QPC electrons and the measured system is dominated by the 
electrostatic potential, the QPC acts as the charge detector. Another example of the 
ballistic mesoscopic detector is the magnetic analog of the QPC based on the ballistic 
motion of the magnetic flux quanta (fluxons) through a one-dimensional channel, the 
role of which is played by the Josephson transmission line (JTL) [T^J. The scattering 
potential Uj{x) for the fluxons in the JTL is created by the magnetic flux or current, 
and the JTL detector can be used for measurements of superconducting flux qubits. In 
the JTL detector, the fluxons can be injected into the JTL individually providing control 
over the individual scattering events. 

The detector model in which the output information is contained in the transport 
current flowing between the two reservoirs applies to many of the mesoscopic detectors 
(see Sec. 4). There are several reasons for this. One is the strong (in the tunnel limit, 
exponential) dependence of the scattering amplitudes on parameters of the scattering 
potential that leads to sufficiently large sensitivity of the detector to the measured system. 
Another, more important, is the fact that the scattering dynamics contains strongly 
divergent transmitted and reflected trajectories that create easily detectable different 
outcomes of measurement. This feature of scattering is not easily reproducible in other 
types of the dynamics 14 . Finally, the transport between large reservoirs makes it 
possible to repeat scattering events at a certain rate amplifying the results of scattering 
of one particle. 

In general, the process of quantum measurement can be understood as creation of 
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Fig. 2. - Measurement dynamics of a ballistic mesoscopic detector. The wavepacket of a particle 
with momentum k is scattered by the potential Uj(x) controlled by the measured system. The 
scattering potential and the transmission/reflection amplitudes tj(k), rj(k) contain information 
about the state \ j) of the system. 



an entangled state of the measured system and the detector as a result of interaction 
between them. The states of the detector are classical and suppress quantum super- 
position of different outcomes of measurement. The two consequences of this process 
are the acquisition of information about the system by the detector, and "back-action" 
dephasing of the measured system - see, e.g., Because of the system-detector 

entanglement, finding a given detector output provides some indication of what state the 
measured system is in. On the other hand, the same entanglement means that quantum 
coherence among the states of the measured system is suppressed. This implies that there 
is a close connection between the information acquisition and back-action dephasing. In 
the optimal situation, the rate W with which the detector obtains information about the 
system and the dephasing rate T are the same. Of course, the detector can always intro- 
duce some parasitic dephasing into the system dynamics, so that in general W < T. In 
view of this inequality, the detector with W = T is called "ideal" or "quantum limited" . 
If the detector is far from being quantum-limited, it destroys quantum coherence in the 
measured system long before it provides information that can be used to select specific 
outcomes of measurement. Because of this, only the detectors that are close to being 
quantum-limited can give rise to non-trivial wavefunction reduction. 

2T. Back-action dephasing rate. - Measurement dynamics with the ballistic meso- 
scopic detector is illustrated in Fig. [21 For the ballistic detector, the detector-system 
entanglement arises as a result of scattering (Fig.[21), and the rates of information acqui- 
sition and back-action dephasing can be expressed in terms of the scattering amplitudes 
[T7|. To do this, we consider evolution of the density matrix p of the measured system in 
scattering of one particle. For simplicity, the Hamiltonian of the system itself is assumed 
to be zero (e.g., A = in the example of Fig. QJ, and the system evolution is caused 
only by the interaction with the detector. The evolution of p is obtained then from the 
time dependence of the total wavefunction of the scattered particle and the stationary 
wavefunction y\ Cj \j) of the measured system: 
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Here ip{x,t = 0) is the initial wavcfunction of the particle injected in the scattering region 
from the reservoir, and its time evolution tpj(x, t) depends on the realization Uj(x) of the 
potential created by the measured system. Tracing over the detector, i.e., the scattering 
wavefunction, one gets from Eq. (Q: 

(2) pij = c t c* — > CiC* J dxipiix, t)ipj(x, t) . 

Qualitatively, the time evolution in Q describes propagation of the initial wavepacket 
towards the scattering potential and its subsequent separation in coordinate space into 
the transmitted and reflected parts that are well-localized on the opposite sides of the 
scattering region. At time t > t sc , where t sc is the characteristic scattering time, the 
separated wavepackets move in the region free from the /-dependent potential and the 
unitarity of the quantum-mechanical evolution of ^pj{x 1 1) implies that the overlap of the 
scattered wavcfunctions in Eq. becomes independent of t. This overlap can be directly 
found in the momentum representation: 

(3) j dx*p i {x,t)ip*(x,t) = Jdkibik^wqw+nWrXk)], 

where b(k) is the probability amplitude for the injected particle to have momentum k in 
the initial state. Equations @ and © show that the diagonal elements of the density 
matrix p do not change in the scattering process: 



(4) J dk\b(k)n ti (k)\ 2 + \ri(W] = l, 

while the off-diagonal elements are suppressed by the factor 



(5) 



J dk\b(k)\ 2 [t t (k)t* (fc) + n {k)r* (k)] 



< 1 . 



The inequality in this relation follows from the Swartz inequality for the scalar product 
in the Hubert space of "vectors" \b(k)\ ■ {tj(k),rj(k)} of the unit length 10}. Suppression 
of the off-diagonal elements of p is manifestation of the back-action dephasing of the 
measured system by the detector. Assuming that the particles arc injected from the 
reservoir with frequency / and combining the suppression factors JSJ for the successive 
scattering events, we obtain the dephasing rate as 



r« = -/In 



dk\b{k)\ 2 [U(k)t*{k) + n(k)r*(k)} 



If the scattering amplitudes do not depend on momentum k in the range of momenta 
limited by |6(fc)| 2 , the back-action dephasing rate becomes independent of the form of 
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initial wavepacket of injected particles: 

(6) r« = -/in|t,t; + r i r;| . 

Equation © is the general expression for the back-action dephasing rate of a ballistic 
mesoscopic detector with momentum- (and energy-) independent scattering amplitudes. 
It is valid, in particular, for the QPC detector at low temperatures T -C eV, if the 
injection frequency / is taken to be equal to the "attempt frequency" eV/h with which 
electrons are incident on the scattering region 24 . 

In the linear-response regime, when the changes in the scattering amplitudes with j 
are small, the limiting form of Eq. © was obtained in ^1 ^| [3U1 1211 |2~2"] . In the tunnel 
limit tj — > 0, Eq. © reduces to 

(7) r^ = (//2)|tw/, 

where 

(8) tj^ltjle^, ^sargfe/rj). 

When the phases 4>j can be neglected, Eq. Q reproduces earlier results (2H| for the 
dephasing by the QPC detector in the tunnel limit. As we will see in Sec. 4, Eq. Q can 
be obtained in this limit under more general assumptions, and describes large number of 
different mesoscopic detectors. 

2'2. Information acquisition rate. - As was discussed above, the back-action dephasing 
is only one part of the measurement process. The other part is acquisition by the detector 
of information about the state of the measured system. In the case of ballistic detector, 
the information is contained in the scattering amplitudes of the incident particles, and the 
rate of its acquisition depends on specific characteristics of the amplitudes recorded by 
the detector. One of the simplest possibilities in this respect, realized, e.g., in the QPC 
detectors, is to record the changes in the transmission probabilities which determine 
the magnitude of the particle current through the scattering region. (Alternatively, 
one could modify the scattering scheme by forcing the scattered particles to interfere, 
and in this way use the phase information |17j in the scattering amplitudes.) The rate 
of information extraction from the current magnitude, i.e., the rate of increase of the 
confidence level in distinguishing different states \j), can be calculated simply by starting 
with the transmission/reflection probabilities Tj — \tj\ 2 and Rj = 1 — Tj when the 
measured system is in the state Since successive scattering event are independent, 
the probability Pj(n) to have n out of N incident particles transmitted, is given by the 
binomial distribution Pj (n) = C^T" '~ n . The task of distinguishing different states \j) 
of the measured system is transformed by the detector into distinguishing the probability 
distributions Pj{n) for different js. Since the number N — ft of scattering attempts 
increases with time t, the distributions Pj{n) become peaked successively more strongly 
around the corresponding average numbers TjN of transmitted particles. The states 
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with different probabilities Tj can be distinguished then with increasing certainty. The 
rate of increase of this certainty can be characterized quantitatively by some measure 
of the overlap of the distributions Pj(n). While in general there are different ways to 
characterize the overlap of different probability distributions |25) . the characteristic which 
is appropriate in the quantum measurement context is closely related to "fidelity" 

in quantum information [23]: J2n\P i ( n )Pj( n )] 1 ^ 2 ■ '^ ne ra * e 01 information acquisition can 
then be defined naturally as [T7|: 

(9) Wy =-(l/*)lnX;[Pi(n)p J -(n)] 1/a . 

n 

Using the binomial distribution in this expression we get: 

(10) Wij = -f\n[{TiTj) 1/2 + {RtRj) 1 ' 2 ] . 

Equation gives the information acquisition rate by ballistic mesoscopic detectors. 
Comparing Eqs. © and (|10|l we see that for this type of the detector, in accordance with 
the general understanding of quantum measurements, the back-action dephasing rate and 
information acquisition rate satisfy the inequality 

(11) wa < r y . 

Equality in this relation gives the condition of the quantum-limited operation of the 
ballistic mesoscopic detector under the assumption of energy-independent scattering am- 
plitudes. It holds if 

(12) = ^ , 

where the phases <fij are defined in Eq. ||SJ). Condition l|12|) has simple interpretation as 
the statement that there is no information on the states \j) in the phases of the scat- 
tering amplitudes. Deviations from Eq. i|12|) mean that the phases contain information 
about the measured system which is lost in the detection scheme sensitive only to the 
transmission probabilities Tj. In this case, the information loss in the detector prevents 
it from being quantum-limited. In practical terms, the simplest way to satisfy condition 
(|12|) is to make the scattering potential symmetric Uj(—x) = Uj(x) for all states \j). 
The unitarity of the scattering matrix for this potential implies then that <j>j = ir/2 for 
any j, and Eq. I|12f> is automatically satisfied. If the detector is quantum- limited, it can 
demonstrate non-trivial wavefunction reduction. 

2'3. Conditional evolution. - Quantitative description of the wavefunction reduction 
due to interaction with a detector can be formulated as "conditional" evolution, in which 
dynamics of the measured system is conditioned on the observation of particular out- 
come of measurement. In the axiomatic approach, wavefunction reduction is formalized 
together with the dynamic evolution as "quantum operation" |27j , arbitrary linear trans- 
formation of the system density matrix satisfying physically motivated axioms. In this 
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approach, a detector is characterized by a set of positive operators which correspond to 
all possible outcomes of measurements with this detector, or "positive operator valued 
measure" (POVM) [2Hj. In practice, for any real specific detector, it is clear what the pos- 
sible classical outcomes of measurements are, and the emphasis is then on development 
of dynamic equations that would describe evolution of the measured system conditioned 
on a given detector output. Since the different outcomes of the detector evolution are 
classically distinguishable, it is meaningful to ask how the measured system evolves for 
a given output. Such a conditional evolution of the measured system describes quantita- 
tively the wavefunction reduction in the measurement process with a particular detector 
(see, e.g., |3!5][j|nilIfflEI|)- In the case of ballistic mesoscopic detectors, each act of particle 
scattering represents an elementary measurement process. Since the particle trajectories 
that correspond to different outcomes of scattering: transmission through or reflection 
from the scattering region are strongly separated, these outcomes should be considered 
as non- interfering classical events. Although the absence of quantum coherence between 
these two outcomes is an assumption of the conditional approach, this assumption is 
very natural. Propagation of the scattered particles in different reservoirs of the detector 
entangles them with different environments, the process that very efficiently suppresses 
their mutual quantum coherence |33| . While this "common-sense" assumption of absence 
of quantum coherence between different outputs of a realistic detector is sometimes con- 
sidered unsatisfactory from an abstract point of view |341 , it can be given a fairly 
rigorous description in terms of decoherence in open quantum systems - see, e.g., |36J. 

Quantitatively, conditional equations are obtained by separating in the total wave- 
function the terms that correspond to a specific classical outcome of measurement and 
renormalizing this part of the wavefunction so that it corresponds to the total probability 
of 1 1161 ITS] . In the ballistic detector, there are two classically different outcomes of 
scattering, transmission and reflection, for each injected particle. This means that the 
wavefunction of the measured system should be conditioned on the observation of either 
transmitted or reflected particle in each elementary cycle of measurement. The evolution 
of the total wavefunction "detector+measured system" as a result of scattering of one 
particle is described by Eq. • Under the assumption of energy- independent scattering 
amplitudes, momentum and coordinate dependence of the states of the scattered par- 
ticles in the detector is the same for different states and can be factored out from 
the total wavefunction. The evolution of the measured system can then be conditioned 
on the transmission/reflection of a particle simply by keeping in Eq. Q the terms that 
correspond to the actual outcome of scattering in the form of the appropriate scattering 
amplitudes. If the particle is transmitted through the scattered region or reflected from 
it in a given measurement cycle, amplitudes Cj for the system to be in the state \ j) change 
then, respectively, as follows: 

(13) Cj -> tjCjf [EjlcjfTj] 1/2 , cj -> TjCjl [Y,j\ C j\ 2 Rj\ 1/2 . 

We see that the expansion coefficients Cj of the system's wavefunction are changing 
in conditional evolution despite the initial assumption that the system Hamiltonian is 
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zero. This is unusual from the point of view of the Schrodingcr equation, and provides 
quantitative expression of reduction of the wavefunction in the measurement process. 

Also, it should be noted that the transformations l|13fl do not decohere the measured 
system despite the back-action dephasing by the detector discussed above. To understand 
this, one should note that as in the case of any dephasing, the back-action dephasing 
can be viewed as the loss of information. For the quantum-limited detection, the overall 
evolution of the detector and the measured system is quantum-coherent and the only 
possible source of the information loss is averaging over the detector evolution. This 
means that the back-action dephasing arises as the result of averaging over different 
measurement outcomes |32) . and specifying definite outcome as done in the conditional 
dynamics removes all losses of information and eliminates the dephasing. 

Equations i|13|) can be applied directly to the detectors which provide control over 
scattering of individual particles, e.g. to the JTL detector [T3J, where it makes sense to 
discuss changes in the wavefunction of the measured system induced by one scattering 
event. In some detectors, however, such a control over individual scattering events is not 
fully possible. For instance, in the QPC detector, the picture of individual scattering 
events leading to Eqs. and (|10|l for the back-action dephasing and information rates is 
strictly speaking valid only on the relatively large time scale t 3> h/eV, when the typical 
number of electron scattering attempts is larger than 1. One can generalize Eq. (|13[1 to 
this situation by considering the time interval t which includes a number N = ft > 1 
of scattering attempts, where for the QPC / = eV/h. Combining the transformations 
(|13|) one can see that observation of any sequence of transmission/reflection events that 
includes n transmissions and N — n reflections changes the wavefunciton as: 



regardless of the specific ordering of these events. This equation includes as a particular 
case Eq. (1131 which follows when N — 1. Since the wavefunction obtained as a result 
of transformation (|14|l is the same for all sequences with the same total number 
n of transmitted particles, one can distinguish all scattering outcomes only by n. For 
each of the N + 1 outcomes with different n the wavefunction is transformed according 
to Eq. I|14|) . This means that the wavefunction reduction has the form (|14|l indepen- 
dently of whether the detector suppresses quantum coherence between all the sequences 
of scattering outcomes or only between the states with different total numbers of trans- 
mitions/reflections. Transformations l|14|) can be used to study quantitatively unusual 
manifestations of the wavefuntion reduction in the mesoscopic solid-state qubits. 

3. Tunneling without tunneling: wavefunction reduction in a mesoscopic 



(14) 




(N-n) 




qubit 



Probably the simplest example of the counter-intuitive features of the wavefuntion 
reduction in mesoscopic qubits arises from a question whether a quantum particle can 
tunnel through a barrier which has vanishing transparency? Immediate answer to this 
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question is "no" as follows from the elementary properties of the Schrodinger equation. 
It seems basically the tautology to say that if the tunneling amplitude is zero (e.g., 
the barrier is infinitely high) the tunneling is suppressed. Of course somewhat more 
careful consideration reminds that evolution according to the Schrodinger equation is 
not the only way for a state of a quantum particle to change in time. Changes in the 
particle state can also be caused by the wavefunction reduction, which, as discussed 
also in the Introduction, can in principle violate any feature of dynamic evolution of a 
quantum system. This can be expressed quantitatively through the "Bell" inequalities 
generalizing the classic Bell's inequalities which quantify violation of the no-action-at- 
a-distance principle. For measurements of a mesoscopic qubit (Fig.^), the peculiarities 
of quantum dynamics of the system originate from the possibility of quantum-coherent 
uncertainty in the position of the charge or flux between the two basis states \j) of the 
qubit. In the regime of coherent oscillations of the qubit (A/0), this uncertainty gives 
rise to several "temporal" Bell inequalities EH1 E01 EH • In this Section, we discuss 
a sequence of quantum transformation that is centered around the qubit dynamics with 
suppressed tunneling, A = 0. The transformations lead to the associated Bell- type 
inequality, which quantify violation of the fundamental intuition of many solid-state 
physicists: charge or magnetic flux can not tunnel through an infinitely large barrier. The 
sequence of transformation includes a measurement done on the qubit and shows that 
the wavefunction reduction can indeed violate this Schrodinger-equation-based intuition, 
and a particle can be transferred through an "impenetrable" barrier in the process of 
quantum measurement. 

The required manipulations of the qubit state are close to those in current experiments 
on coherent oscillations and more complex dynamics of mesoscopic qubits. Although the 
discussion in this Section applies in general to all types of qubits, the physics content of 
the wavefunction reduction is more striking for the semiconductor 111! I12| or super- 
conductor 13211131311131101 charge qubits, or for the flux qubits g7J EHl ESI In this 
case, the basic set-up is equivalent to the one shown in Fig. ^ The two basis states \j), 
j = 1,2, of the qubit differ by some amount of magnetic flux or by an individual charge 
(electron charge e in semiconductor quantum dot qubit, or Cooper-pair charge 2e in a 
superconductor qubit) localized on the opposite sides of a tunnel barrier. The states are 
coupled by the tunnel amplitude A > 0. At the point of resonance, when the bias energy 
e between the two basis states vanishes, the qubit Hamiltonian reduces to 

(15) H = -Aa x , 

and describes quantum coherent oscillations with frequency 2A. The tunnel amplitude 
A and the bias energy e are assumed to be controlled externally. 

The main element of the sequence of qubit transformations consists in preparing 
a superposition (for simplicity, symmetric superposition, a x — 1) of the qubit basis 
states, then switching off the tunnel amplitude A and performing quantum-limited but 
weak measurement of the qubit position in the a z basis. The idea is to prove that the 
measurement-induced transfer of the qubit wavefunction between the two qubit states at 
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A = gives not only the changes of the probability (representing our knowledge about 
the qubit position) but the actual transfer of charge or flux through the infinitely large 
barrier. In order to do this, we include the measurement-induced wavefunction transfer 
as a part of the cyclic transformation, the other part of which is known to transfer 
charge or flux. In the ideal situation, when the operations are precise and there is no 
intrinsic decoherence, the cycle should lead to precisely the same initial state cr x = 1, 
making it possible to conclude that the measurement part of it transferred the qubit 
state through the barrier with vanishing tunneling amplitude. In the presence of external 
perturbations, there will be a probability p( — ) to find the qubit not in the initial state. 
One can, however, derive a condition in the form of inequality on p(~^ , violation of which 
shows that this observation cannot be explained within the assumption of some initial 
classical probability distribution over the qubit basis states. This means that explanation 
of the observed violation should necessarily involve the transfer of the qubit state through 
suppressed barrier by the process of the wavefunction reduction. 

In detail, the starting point of the sequence of transformations is the ground state of 
the Hamiltonian l|15|) 

(16) |^) = (|1) + |2))/V2, 

in which the qubit will find itself because of the unavoidable, but assumed to be weak, 
relaxation processes, if the Hamiltonian (|15|l is kept stationary for some time. Starting 
from this state, the tunnel amplitude A is abruptly switched off, A = 0. The rate of 
this process is not important in the case of the Hamiltonian (|15|l . since the final state 
will coincide with regardless of how slowly or quickly A is switched off. However, in 
the presence of some parasitic residual bias e, the rate of variation of A should be much 
larger than e to preserve the state (11611 at the end of the switching process. Next, as the 
first step of actual transformations of the qubit state l|16|) . the weak quantum- limited 
measurement of the a z operator is performed on this state. The result of this operation 
does not depend on the specific model of measurement, as long as it is quantum-limited. 
For such a measurement, we know that specifying the detector output n leaves the qubit 
in a pure state which is obtained from by increased degree of localization in the a z 
basis because of the information about a z provided by the measurement: 



(17) 



|V>i)=a„|l)+/3„|2) 



As a suitable model of this measurement one can use the measurement with a ballistic 
detector discussed above. In this case, the detector output is the number n of transmitted 
particles in N = ft scattering attempts. In this case, according to Eq. I|14|l 



f nJN-n) 



T n R (N-n) 



n r>(N-n) 



T?R 2 



1/2 



Pn = ~ 



f nJN-n) 
l 2'2 



T n R {N-n) 



njy(N-n) 



T?R 2 



1/2 



Using the condition i|12|) of the detector ideality one can see that the relative phase £ of the 
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Fig. 3. - Probability amplitude a n l|18fl of finding the qubit in the state |1) as a function of the 
observed number n of particles transmitted in TV = 10 attempts through the ballistic detector 
with transmission probabilities 21 = 0.8 and T2 = 0.4. The detector measures the state 116H . 



coefficients a n and (i n is independent of the detector output n: £ = iV[arg(ti) — arg(i 2 )]- 
It can be viewed then as renormalization of the qubit bias e due to the detector-qubit 
coupling, and can be compensated for by the bias shift Se = /[arg(ii) — arg^)] during 
the period of measurement. The coefficients a, (3 have then the following form: 



(18) 



where 





1/2 
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can be interpreted as the relative probability for the qubit to be in the state \j) for a 
given detector outcome n. In the situation when the detector provides no information 
on <r z , e.g. if T\ = T2, the coefficients (|18|) are unchanged from their initial values (|16fl . 
Otherwise, the probability amplitude is shifted in the direction of the more probable 
state: the amplitude (|18l) of one qubit state is increased in comparison with l|l(j|) if the 
observed n is closer to the value rij — TjN characteristic for this state than the other 
state. As an illustration, Fig. [3] shows the amplitude a n for TV = 10, T\ — 0.8, and 
T-2 = 0.4, as a function of n. One can see that a n maintains its original value l/y/2 from 
Eq. (|16fl for n's roughly in the middle between the two characteristic values n\ — 8 and 
712 = 4. For n smaller than or close to n2, a n decreases from its original value, for n close 
to or larger than n±, a n approaches 1. Such a shift due to the wavefunction reduction is 
the central part of the transformation cycle. 

The remaining steps of the cycle aim at returning the qubit to its initial state (|16f) . 
To do this, one needs to transfer back the charge or flux that was transferred in the 
wavefunction reduction process leading to the state l|18|l . This is achieved by creating for 
some time the non-vanishing tunneling amplitude, i.e. realizing a fraction of a period of 
the regular coherent oscillations in which the charge or flux goes back-and-forth between 
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the two qubit basis states. In the most direct way, this can be done if the qubit structure 
makes it possible to create non- vanishing phase of the tunnel amplitude A'(£) (e.g., in the 
superconducting qubits, where the tunnel amplitude can be controlled through quantum 
interference, producing any complex value of this amplitude). In this case, the state (|18|) 
can be returned back directly into the initial form (|16(l if argA' = ir/2. In the diagram 
(Fig. ^Ji) in which the qubit states are represented in the language of spin-1/2, i.e. 

|^} = cos(0 n /2)|l}+sm(0 n /2)|2>, 

such a tunneling amplitude corresponds to rotation around the y axis. The diagram in 
Fig-EJi shows then directly that the rotation around y axis turning into | V^o) should 
have the magnitude: 



(19) / \A'(t)\dt/h= (%/2-6 n )/2, 
where 

6 n = 2tan- 1 (/3„/a„) = 2t.au- 1 [(T 2 /T 1 ) n / 2 (R 2 /Ri^ N - n)/2 

If the qubit structure allows only for the real tunnel amplitude A (the situation 
that can be expected, e.g., in semiconductor quantum dot qubits), the y-axis rotation 
R y = cxp{—ia y j \A'(t)\dt/h} ()19(1 can be simulated in three steps in which the rotation 
R x = e~x.p{—ia x f A(t)dt/h} around the x axis of the same magnitude l|19fl is preceded 
and followed by the rotations around the z-axis: 

(20) R y = U- X R X U , U = exp{icr 2 7r/4} . 

The z-axis rotations can be created by the pulses of the qubit bias: J e(t)dt/h — ±ir/4. 
The three-step sequence (|2U|1 can be simplified into two steps (Fig. 0Jd) by changing 
the order of rotations: first, the x-axis rotation by 7r/4 (opening tunneling A(t) for 
appropriate interval of time) followed by one z-axis rotation: 



(21) 



j A(t)dt/h = 71-/4, J e(t)dt/h = (tt/2 - 9 n )/2 . 



All these versions of the transformation cycle bring the qubit to its initial state (|16fl • In 
all cases, completion of the cycle that started with a shift of the wavefunction amplitudes 
due to the state reduction involves a part of a period of coherent oscillations which 
reverses this shift. Coherent qubit oscillations are known to actually transfer the charge 
of flux between the two qubit states. Since the cycle as a whole is closed, this fact 
shows that the changes in the qubit state caused by the wavefunction reduction can not 
be interpreted only as the changes in our knowledge of probabilities of the state of the 
qubit, but involve the actual transfer of the charge or flux in the absence of the tunneling 
amplitude. 



14 



D.V. Averin 



(a) 



z 



(b) 



z 




1^1 > 




IV x 



y 



Fig. 4. - Diagram of the two possible transformations of the qubit state in the spin representation 
returning the state \ipi) to \ipo) after the measurement-induced state reduction \ifio) — ► l^i)- 
(a) Direct one-step i/-axis rotation l|19|l . (b) Projection on the z — y plane of the two-step 
transformation 12 II with the same end result. 



To see this more quantitatively, one can derive the Bell-type inequality, violation of 
which should show that understanding of the state reduction solely in terms of probability 
changes can not be correct. The inequality is obtained by assuming that the process of 
switching off the tunneling amplitude A in the beginning of the transformation cycle 
does not leave the state (tTtjll unchanged but instead localizes the qubit in one of the 
basis states. This means that the process produces an incoherent mixture of the qubit 
states with some, in general unspecified, probability p to be in the state |1). This 
process would provide then an alternative, classical description of the evolution during 
the measurement process. In this description, the qubit state is "objectively" well defined, 
but is unknown to us, and the measurement gradually provides information about this 
unknown state. The measurement would only change the probabilities we ascribe to the 
two qubit states, but not the state itself, and in particular would not transfer the charge 
of flux. One should then see how well this classical description can mimic the quantum 
result of the transformation cycle described above. A convenient way of making this 
comparison is provided by the probability of ending the cycle in the wrong state. The 
unperturbed quantum evolution should end up in the initial state u x = 1, whereas the 
same transformation cycle performed on the classical initial state will always have a finite 
probability p^' of ending in the state a x = —1. 

This probability is found by applying the transformations not to the state <|16f) but 
to the incoherent state with the density matrix 



Po=p|l)(l| + (l-p)|2)(2|. 



In this case, the measurement changes only the probability p in this expression. Similarly 
to Eq. ljT%|) . if the detector gives the output n, the density matrix of the system is 
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All versions l|19fl - (|21|l of the transformation cycle produce the same probability p( - ) of 
being in the state a x = — 1, when applied to the density matrix p\ (|22J) . For instance, 
one can see directly that in the density matrix R~ l p\R y obtained from p\ by rotation 
(|T§|) . the probability p ( ~) is: 

p(-) = !™ . 

(wi + w 2 )(pwi + (1 - p)w 2 ) 

Minimizing this expression with respect to p, we see that the minimum probability of 
finding the state a x = — 1 in the classical case is 

(23) p(-> = 1 - 2J . 

Wl + W 2 

Instead of looking for the minimum with respect to p, one can adopt a natural additional 
assumption that when the tunneling amplitude is switched off, the qubit localization 
process can only be symmetric, since there is no reason to prefer one qubit state to 
another. In this case, p = 1/2, and we obtain somewhat different expression for the 
probability p^ with qualitatively similar properties: 



(24) p 



(wi + w 2 ) 2 



This expression would also be obtained if the qubit wavefunction would be reduced to 
the density matrix p\ (|22|l during the measurement, not in the process of suppression of 
the tunneling amplitude. 

Equations (|23[) and l|24l) show that in order to distinguish the quantum coherent 
evolution (for which p(~) = 0) and incoherent evolution with non-vanishing probability 
p(~\ it is important to employ weak measurement. If the measurement is projective, 
i.e. if one of the probabilities Wj is zero so that the measurement completely reduces the 
qubit state to one of the basis states, then p^ — and it is impossible to distinguish 
the two types of evolution. This conclusion should be independent of the specific form 
of the employed transformation cycle, since projective measurement is always expected 
to fully separate different components of the initial state of the measured system and 
completely suppress quantum coherence between them. 

The discussion above means that observation of the probability of the state <j x = — 1 
smaller than p^', 

(25) p(a x = -1) < 

at the end of the transformation cycle proves that all transformations in this cycle, 
including the wavefunction reduction, are quantum coherent. Combined with the non- 
vanishing transfer of charge or flux during the "oscillation" step [ ()19f) - l|21|l ] of the 
cycle, this fact implies that the wavefunction reduction induces similar transfer across 
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the tunnel barrier separating the qubit basis states even if the corresponding tunnel 
amplitude is zero. 

It is important to note that this counter-intuitive feature of the wavefunction re- 
duction does not contradict the fact that all dynamic properties of the measurement, 
including the back-action dephasing and information rates can be calculated from the 
dynamic detector model without any reference to the wavefunction reduction. While 
the dynamic properties are average characteristics of the detector, the wavefunction re- 
duction appears if one considers separately individual outcomes of measurement. In the 
transformations discussed above, this separation is achieved by introducing the feedback, 
operations on the measured system which depend on the specific measurement outcome. 
In ballistic detectors, the measurement outcomes are distinguished by the number n of 
transmitted particles, and to see the wavefunction reduction one needs to distinguish 
individual particles, and is done, e.g., in the electron counting experiments [5U|. For the 
detectors for which distinguishing individual transmitted particles can be problematic 
(e.g., the QPC detector), the allowed uncertainty in n should be smaller than the width 
5n of the transition region in the n-dependence of the wavefunction amplitudes of the 
measured qubit - see Fig. [21 Since the width of this region can be estimated roughly 
as Sn ~ 1/ \n{TiR2 / T2R1) , the uncertainty in n can be compensated for by making the 
difference between transmission probabilities Tj smaller, thus increasing 8n. The limit to 
this increase is set by the decoherence processes in the measured system which make it 
impossible to increase the measurement time of the detector beyond the coherence time 
of the system without losing the non-trivial character of the measurement dynamics. 

4. Tunneling detectors 

So far, the discussion was based on the ballistic model of the mesoscopic detector, 
in which the measured system controls ballistic motion of some particles between the 
two reservoirs (Fig. |2J. If one assumes that the particle transmission probabilities are 
small, Tj <C 1, and the transfer processes between the reservoirs can be described in the 
tunneling approximation, specific nature of the detector transport becomes irrelevant. In 
this case, the range of applicability of the detector model can be extended significantly to 
include the detectors in which it is not possible to identify regions of ballistic transport, 
but which are still based on the very similar dynamic principle: control by the measured 
system of transport between the two reservoirs. Examples of such "tunneling" detectors 
include the superconducting SET electrometer |511 1521 , normal SET electrometer in 
the co-tunneling regime [231, or dc SQUID magnetometer (see, e.g., (^ESJISniEl) used 
for measurements of superconducting qubits. The aim of this Section is to show briefly 
that the measurement properties of this type of mesoscopic detectors coincide in essence 
with those of the ballistic detectors. 

Since the measured system controls the tunneling amplitude t of particles in the 
detector, this amplitude should be treated as a non-trivial operator acting on the mea- 
sured system. The detector tunneling can be described then with the standard tunnel 
Hamiltonian, the transfer terms in which are split into a product of operators of the 
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measured system and the detector. In this case, the tunnel Hamiltonian describes the 
detector-system coupling and can be written as 



where are the operators that describe the detector part of the tunneling dynamics, 
e.g., creation of excitations in the detector reservoirs when a particle tunnels, respec- 
tively, forward and backward between them. Inclusion of the operators £, £t means, 
therefore, that the Hamiltonian 126(1 makes it possible to describe the detectors in which 
the tunneling processes are strongly inelastic. This fact, however, does not prevent cor- 
rect account of elastic transport in the case of ballistic detectors. Qualitative reason 
for this can be seen easily using as an example the QPC detector. While scattering of 
individual electrons in the QPC is elastic, in the tunnel limit, electron transfer between 
the two electrodes of the QPC can also be viewed as creation of electron-hole excitation 
in the electrodes, with an electron removed from a state below the Fermi level in one 
electrode, and transferred to a state above the Fermi level in the other. Accordingly, as 
we will see later, Hamiltonian (|26[) leads to the evolution equations that coincide in the 
tunnel limit with those obtained above in the ballistic case. 

Under the assumption that the detector tunneling is weak, the precise form of the 
internal detector Hamiltonian is not important and dynamics of measurement is defined 
by the correlators of the operators 



Here the angled brackets denote averaging over the detector reservoirs which are taken 
to be in a stationary state with some fixed number of particles in them and the density 
matrix p^: (...) = Tr d{---Pd} ■ The correlators (|27|) set the scale T± = 2Re7± of the 
forward and backward tunneling rates in the detector. 

A reasonable tunneling detector should satisfy some additional assumptions related to 
the fact that its output should be classical in order to provide a complete measurement 
dynamics. Similarly to the ballistic detector, the output information in the tunneling 
case is contained in the number n of the particles transmitted between the detector 
reservoirs. For this number to behave classically, the correlators (£(£)£), (^{t)^) that 
do not conserve the number of tunneling particles should vanish. Another consequence 
of the assumption of classical detector output is that the energy bias AE for tunneling 
through the detector should be much larger than the typical energies of the measured 
system. In this case, one can neglect quantum fluctuations of the detector current in the 
relevant frequency range that corresponds to frequencies of evolution of the measured 
system. If one assumes in addition that all other characteristic frequencies of the detector 
tunneling are also much larger than those of the measured system, the functions £(i), 
£' (i) in Eq. I|27|) are effectively ^-correlated on the time scale of the system. 

Vanishing correlation time in the correlators l|27|l makes it possible to write down 
simple evolution equations for the density matrix p of the measured system. To describe 



(26) 



H T = % + tV , 



(27) 
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the system dynamics conditioned on particular outcome of measurement, we also keep in 
the evolution equation the number n of particles transferred through the detector. Since 
the correlators that do not conserve n vanish, only the terms diagonal in n are important. 
In the interaction representation with respect to the tunnel Hamiltonian H26[l . the density 
matrix p{i) is given by the standard expression: 

(28) p(t) = (Spp D S^), S = Texp{-if dt'H T {t')}. 

For ^-correlated operators in Eq. I|27|). one can see that the full perturbation expansion 
of Eq. H28J1 in Ht is equivalent to the evolution equation for p(t) that follows from the 
lowest-order perturbation theory. Keeping track of the number n of particles transferred 
through the detector, we get: 

(29) p {n) = r+ ty n_1) * + r_ i P (n+1) P - { 1+ W + 7_t t tV n) - p {n) (7+^ + 7* tH) . 

Since the tunneling amplitude t is a function of some observable of the measured 
system, there is a system of eigenstates |j) common to the operators t and v: 

i\j)=tj\j), $\j) = t*\j), 

where tj is the detector tunneling amplitude when the measured system is in the state 
\j). It is convenient to write the evolution equation (|29|l in the basis of states |j): 

(30) 4 n) = -(i/2)(r + + r_)(|£ 2 | 2 + M 2 )p? } + r_^t>J +1) 

+T + t*t 3 p^- 1) -i[5H, p(%j, 

where 

(31) ^ = Im( 7++7 _)|t J | 2 |j)(j| 

is the renormalization of the Hamiltonian of the measured system due to its coupling to 
the detector. 

If one omits the term 5H which can be combined with the internal Hamiltonian of the 
measured system, Eq. (|30ll can be solved in n by noticing that it coincides in essence with 
a recurrence relations for the modified Bessel functions /„ |58|. In order to interpret n as 
the number of particles transferred through the detector during the time interval t, we 
solve this equation with the initial condition Pij(t = 0) = Pij(0)S n ,Q. The corresponding 
solution is: 

o {n) (t) /r \ n/2 r +r 

(32) ^-^y = (^J I n (2t\t l t J \^TTrZ)cM-^-(\t^\ 2 + \tJ\ 2 )t-^n^ J }, 
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where ipij = axg(tit*). As discussed above, the qubit density matrix conditioned on the 
particular measurement outcome n follows from Eq. (|32|) if one selects in this equation the 
terms with given n and normalizes the resulting reduced density matrix back to 1. If one 
of the tunneling rates T + or T_ vanishes, Eq. I|32l) reduces to the usual Poisson distribu- 
tion characteristic for tunneling in one direction. Conditional description of measurement 
in this case was developed in JSj. When both rates are non- vanishing, specifying the 
total number n of the transferred particles does not specify uniquely evolution of the 
detector, since the same n results from the balance between different numbers of parti- 
cles transferred forward and backward. This means that some information is lost in this 
regime and the detector is not quantum-limited (see the discussion below). In contrast 
to the situation with the quantum-limited detectors considered in the previous Sections, 
conditional dynamics that follows from Eq. (|32|) in this case does not preserve the 
purity of the quantum state of the measured system. 

Evolution of the density matrix p averaged over the measurement outcomes n can be 
obtained by either disregarding simply the index n in Eq. I|30() . or directly taking the 
sum over n in (|32[l with the help of a summation formula j58j for the Bessel functions. 
In both cases, equation for the measurement-induced evolution of p is: 

(33) pij = -Tijpij - i(T + - r_)|i;ty| -in ,;,,/',, , 
where 

(34) r ij = (i/2)(r + + r_)\t i -t j \ 2 

is the back-action dephasing rate by the tunneling detector. The last term in Eq. H33fl can 
be viewed as another contribution to renormalization of the energy difference between 
the states |i) and \j), although in general it can not be reduced to the energy shifts of 
individual states [in contrast to the renormalization term (I31[) ] . 

The back-action dephasing rate (|34|l coincides with that of the ballistic detectors 
in the tunneling limit given by Eq. J7J, if we take into account that our discussion of 
ballistic detectors assumed for simplicity that particles are incident on the scattering 
region from only one electrode. It can be seen directly that the difference in the phases 
[see Eq. ©] of the tunneling amplitudes in the two expressions is not essential. The 
reason for this difference is that the tunnel Hamiltonian (|26|) describes explicitly only the 
effects associated with the actual transfer of particles across the detector. It is assumed 
that the other effects of the detector-system coupling, e.g., renormalization of the system 
energy due to reflection processes in the detector, are already accounted for. In terms of 
the scattering amplitudes, this implies that the reflection amplitudes r and f for particles 
incident on the tunnel barrier from the two electrodes of the detector, should satisfy the 
conditions arg(r^r*) = and arg(f#|) = 0. In this case, Eq. Q coincides with l|34|) 
even with account of extra phases ijSJl • This means that the model of the tunnel detector 
considered in this Section is equivalent to that of the ballistic detector in the appropriate 
small-transparency limit. 
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Fig. 5. - The information acquisition rate Wij of the tunneling detector normalized to the 
dephasing rate I\j II34H (for ipij — 0) as a function of time t for several ratios of the forward 
and backward tunneling rates. The time t is normalized to the typical forward tunneling rate 
7 = r+ (T, + Tj)/2. The curves are plotted for the detector transparencies T, = 0.2, Tj = 0.4. 
The dashed lines show the corresponding asymptotic values (I35H . For r+/r_ = 100, the dashed 
line overlaps with the main curve. 



To calculate the information rate Wij of the tunneling detector in the situation when 
the rates of both forward and backward tunneling are non-vanishing (in contrast to 
scattering from one direction discussed for the ballistic detectors), one needs to use in 
Eq. © the diagonal part of Eq. (jMl'li which gives the probabilities Pj(n) — p^' for n 
particles to tunnel when the system is in the state An example of the rate Wy 
defined in this way is shown in Fig. [SJ This figure shows that in general Wij is time- 
dependent and approaches constant value only after a transition period. If the tunneling 
probabilities Tj — \tj\ 2 do not differ very strongly, this transition period is shorter than 
the back-action dephasing time Iy. 

The constant asymptotic values of the information rate can be obtained from the 
asymptotic behavior of the Bessel functions I n (z). Using the standard integral represen- 
tation for I n (z) in Eqs. © and l|32|) . and making use of the fact that the constant rates 
Wij arc determined by the exponential behavior of the integrals at large time t, we find 
directly 

(35) w^ = (r+ + r_)(T, + Tj)/2 - [{t? + r?)r+r_ + lyr^r^ + r 2 _)] 1/2 . 

If the particles tunnel only in one direction, e.g. T_ = 0, Eq. I|35l) reduces to the 
previously known result E2 Wij = T + (^/Ti — y /7 Fj) 2 /2, in which the information 
rate and the back-action dephasing rates coincide, when the phases of the tunneling 
amplitudes satisfy the appropriate ideality condition ifij = 0. In the other limit of small 
difference 2AT between the transmission probabilities Tij = T ± AT, Eq. 1|35|) reduces 
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to 



(36) Wa 



(Ar) 2 (r H 

2T(r+ 



This equation agrees with the general theory of linear measurements (see, e.g., j!5]^ . in 
which it can be interpreted as the rate with which one can distinguish the difference 
AT(r + — r_) between the two detector currents in the presence of the current noise 
with spectral density T(T + + 

In the most typical situation, the two rates r± are both non-vanishing because of the 
finite temperature of the detector electrodes. The electrodes can be in equilibrium 
even when a non-vanishing current is driven between them by finite energy difference 
AE created for the tunneling particles. In this case, the tunneling rates r± are related 
by the detailed balance relation and can be written as r± = To exp(±Ai?/29), where 
To is the typical tunneling rate which can also depend on temperature and energy bias. 
The information rate 1|35[1 then is 



(37) W i:j /To = (Ti + Tj) cosh(AE/2G) - [T? + Tf + 2TiT cosh(A£:/e)] 1 



Comparison of Eqs. H37f) or (|35|) with Eq. 1341) for the back-action dephasing rate (see 
also Fig. [SJ) shows that the detector with temperature ~ AE which creates non- 
vanishing rates F± , is not quantum-limited, Wij < , even if the phases of the tunneling 
amplitudes satisfy the ideality condition = 0. Equation l|37|) can be used to establish 
quantitative condition on the detector temperature necessary for the desired degree of 
the detector ideality. 



5. Conclusion 

Two general models of realistic mesoscopic solid-state detectors have been described 
in this paper. The detectors are based on ability of the measured system to control 
transport current between two particle reservoirs. The models enable detailed analysis 
of the dynamics of the measurement process. Wavefunction reduction is introduced in 
this dynamics through the assumption of suppressed quantum coherence between the 
particle states in different reservoirs. This procedure is very natural and can be justified 
qualitatively within the general approach to decoherence in quantum systems. The main 
element of the justification is the increased level of difficulty of maintaining quantum 
coherence between the states of progressively more complex systems. This fact makes 
the boundary between the quantum and classical domains not very sharp and dependent 
on details of specific measurement set-up. This leads to an interesting question whether it 
is possible to formulate more general and self-consistent conditions defining the boundary 
between the quantum and classical behaviors of dynamic systems. Mesoscopic solid-sate 
structures provide a convenient setting for further studies of this question. 
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* * * 

Part of the discussion in this paper is based on the work done in collaboration with A. 
di Lorenzo, A.N. Korotkov, K. Rabenstein, R. Ruskov, V.K. Semenov, E.V. Sukhorukov, 
and W. Mao. The author would like to thank them, and also G. Benenti, R. Fazio, J.W. 
Lee, F. Plastina, and D.L. Shepelyansky, for collaboration and discusssions of quantum 
measurements. This work was supported by the NSF under grant # DMR-0325551. 
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